Remarks on Flory theory of a self-avoiding chain under cylindrical confinement 
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Much progress in understanding chain molecules often 
thrives on simplification, due to their intrinsic complex- 
ity. Nevertheless, it is hard to exaggerate the success 
and the impact of Flory's brilliant scheme for comput- 
ing the equilibrium size of a swollen polymer chain in a 
good solvent Consider such a swollen polymer car- 
rying N monomers with its end-to-end distance R. The 
trial Flory free energy in d spatial dimensions is then 
expressed as [2] 
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where a the monomer size, f3 — 1/ksT with ks the Boltz- 
mann constant and T the absolute temperature. (Here 
and below, we only consider an athermal case, i.e., the 
excluded volume of each monomer is a d , independent 
of T.) The first term describes chain elasticity of en- 
tropic spring, while the second term represents the mean- 
field energy of (two-body) interaction between monomers 
along the chain [l], Q . The equilibrium chain size or the 
Flory radius, Rp, is obtained by minimizing J- with re- 
spect to R: For d < 4, R F ~ aN v with v = 3/(2 + d) 
(e.g., v — 3/5 for d = 3). This exponent v is rightly 
designated as the Flory exponent. 

As de Gennes correctly pointed out, however, Flory's 
theory benefits from the remarkable cancellation of errors 
(overestimates) of both terms in Eq. [I] 3]. Note that, at 
R = Rp, the Flory free energy scales as 



(2) 



This result is equivalent to stating that the number of 
monomer contacts scales as N 2 /R F ~ N 2 ~ vd (~ N't for 
d = 3 [4|). For d — 3, this is an overestimate in view of 
the more elaborate analysis in Ref . [H, Q . On the other- 
hand, /3^ci=4(Rf) ~ 1 is expected to be asymptotically 
valid. 

Despite the aforementioned limitations, mainly due to 
its simplicity, Flory-type approaches have been extended 
to many other important cases, e.g., linear chains with 
stiffness 0, H| and polymers of various topology in a con- 
fined space (for a review, see Ref. [9]). In particular, the 
widely-used free energy of a linear chain in a cylindrical 
pore has the following form [1(3] 
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where R\\ is the trial chain size in the longitudinal direc- 
tion and D the width of the cylinder. 

Although Eq. [3] produces the correct scaling for the 
equilibrium chain size % - Na{a/D)i its ap- 
plicability beyond the computation of i?|| is severely 
limited by the following two important errors: (i) For 
a chain in equilibrium (R\\ — R\\q), Eq. [3] predicts 
(3 T ~ N(a/D)i, but the correct free energy should 
scale as (3T - N(a/D)i 0, [T3|. (ii) Eq. M results in 
an incorrect effective "Hookian" spring constant k e g of 
the chain, f3k e s ~ 
[3k eS - N^D-ia 
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but the correct scaling is 
|l2j. The two errors are not 
unrelated, and, in fact, arise from the same source. For 
instance, the entropy term of Eq. [3] fails to correctly 
take into account the additional length scale D associ- 
ated with the cylinder - it erroneously assumes that the 
elasticity is not affected by the presence of confinement. 

Here, we propose the following "renormalized" free en- 
ergy for a polymer under cylindrical confinement. 
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where g is the number of monomers inside a blob of di- 
ameter D, i.e., g ~ (D/a)3. Our basic idea is to con- 
sider the confined space inside a cylinder as an effec- 
tive one-dimensional space {d = 1) and, hence, to intro- 
duce the lengthscale D accordingly in Eq. Q] by rescaling 
a -> a' = D and N -> N' = N/g. Then, the first term 
can be understood as the chain being made of N/g sub- 
units ("blobs") of size D, while the second term describes 
the mutual exclusion between neighboring blobs. In other 
words, Eq. [4] represents a self-avoiding chain confined in 
a cylinder of diameter D as a one-dimensional Rouse (or 
Gaussian) chain with an effective step length D, since 
long-distance interactions beyond D (along the contour) 
are suppressed; thus, it is also expected to exhibit the 
Rouse dynamics 11 1. 

Indeed, the above free energy produces not only 
the expected equilibrium chain size Ruq ~ Na(a/ D)~s , 
but also the correct blob-overlapping free energy of 
/3-^conf ~ N/g ~ N(a/D)z (namely, the total number 
of blobs) @, llll . Note that this is identical to con- 
finement free energy. The way neighboring blobs inter- 
act each other is not different from the way they are 
constrained by the confining wall. Importantly, we ob- 
tain the correct effective Hookian spring constant of the 
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chain: (3k e g 



N~ 



~(a/D)3 [11|, [12|. Moreover, using the stretch- 



release argument, the global (slowest) relaxation time of 
the confined chain (in the absence of hydrodynamic ef- 
fects) is then reciprocally obtained as tr ~ N/k e g ~ 
N 2 a» D 2 ~» ~ N 2 a 2 (D/a)3 , consistent with the earlier 
scaling result 11, T3|. 

A natural extension of Eq.[2]to the slit case [l3| would 
be obtained from Eq. [1] with the same rescaling of a — > 
a' = D and N — » iV' = -/V/g, assuming d — 2, as follows 
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where I? is the distance between the two parallel slits 
and g ~ (D/a) 5 ^ 3 . Indeed, this free energy predicts the 
correct equilibrium size of the chain in a slit, i?|| ~ 
aiV 3 / 4 (a/D) 1 / 4 , in agreement with the results of other 
approaches [H, @|- However, its equilibrium free energy, 
/3^ s iit(i?|| = R\\o,D), scales as iV 1 / 2 (a/D)i. Note that 
this is different from the free energy of slit confinement 
/S-^conf ~ AT/ g ~ Nja/D) 3 , which is identical to that of a 
cylinder (see Ref. [3, [l3| and references therein) . The slit 
free energy in Eq. [5] describes the arrangement of blobs 
in a slit, not confinement. 

Why does the Flory approach work better for the cylin- 
der case than for higher dimensions? To see the differ- 
ence, consider the volume fraction (a) of monomers in- 
side a space explored by a self-avoiding chain in equilib- 
rium. For cylindrical confinement, this is given by a = 
Na 3 /D 2 Rp ~ (D/a)~ 4 / 3 , independent of N, whereas, 
both for a slit and for a dilute bulk solution, a — > 
as N — ► oo. It is this special nature of compactness 
of one-dimensional space that explains the extensiveness 
and thus the success of the rescaled Flory approach in 
Eq. Hill- 
Extension of our analysis to semi-flexible polymers and 
chains with non-trivial topology (e.g., branched) is im- 
portant for understanding many recent experiments on 
DNA in nano- and micro-channels [l5| and have non- 
trivial implications for practical applications such as fil- 
tration [16|] . 
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